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We present a theoretical study of the dynamical behavior of a gas made of ultracold fermionic 
atoms, which during their motions can collide with a much smaller number of thermal bosonic 
impurities. The atoms are confined inside harmonic traps and the interactions between the two 
species are treated as due to s-wave scattering with a negative scattering length modeling the '"'K- 
*^Rb fermion-boson system. We set the fermions into motion by giving a small shift to their trap 
center and examine two alternative types of initial conditions, referring to (i) a close-to-equilibrium 
situation in which the two species are at the same temperature (well below the Fermi temperature 
and well above the Bose-Einstein condensation temperature); and (ii) a far-from-equilibrium case in 
which the impurities are given a Boltzmann distribution of momenta while the fermions are at very 
low temperatures. The dynamics of the gas is evaluated by the numerical solution of the Vlasov- 
Landau equations for the one-body distribution functions, supported by some analytical results on 
the coUisional properties of a fermion gas. We find that the trapped gaseous mixture is close to the 
coUisionless regime for values of the parameters corresponding to current experiments, but can be 
driven towards a coUisional regime even without increasing the strength of the interactions, either 
by going over to heavier impurity masses or by matching the width of the momentum distribution 
of the impurities to the Fermi momentum of the fermion gas. 

PACS numbers: 03.75.Ss, 02.70.Ns 



I. INTRODUCTION 

Trapped spin-polarized Fermi gases can be considered 
as non-interacting at ultralow temperatures since s-wave 
collisions are forbidden by the Pauli principle and higher- 
wave collisions are negligible The coUisionality of 
such gases has been increased by mixing it either with a 
fermion gas in a different spin state H, ^ lEi 1^ or with 
a Bose-Einstein condensed gas of bosonic atoms in num- 
bers exceeding those of the fermions by a few orders of 

mag nitude laaJMES- 

However, the coUision rate of 
a fermion-fermion mixture at low temperatures is limited 
by the Pauli blocking of collisions as a result of the occu- 
pation of final states around the Fermi level and, in the 
case of a boson-fermion mixture, the coUisionality can be 
strongly diminished by superfluidity of the condensate 
[T^ . On the other hand, it has been shown by Ferlaino 
et al. jl^ that the number of collisions increases if the 
bosons are thermal and this can be realized by dimin- 
ishing the condensation temperature as for instance can 
be achieved by lowering the number of bosons. It was 
previously pointed out by Amoruso et al. |Q| that the 
presence of a small number of bosonic impurities could 
drastically increase the coUisionality of a Fermi gas, to 
the point of driving it from the coUisionless to the coUi- 
sional regime. However, the Fermi sphere is almost fully 
occupied in ultracold gases close to equilibrium as are 
realized in actual experiments, and the Pauli principle 
will in this case still limit the scattering against impuri- 
ties. The dependence of coUisionality on the presence of 
impurities thus needs detailed investigation. 

In this work we study the coUisional properties of an 
atomic Fermi gas interacting with a few atomic impuri- 



ties. The presence of the impurities induces a damping of 
oscillatory motions of the gas and a shift of its natural os- 
cillation frequencies. We focus on the collision and damp- 
ing rates as functions of the concentration of impurities 
and of temperature for system parameters corresponding 
to *°K-^^Rb fermion-boson mixtures relevant to current 
experiments [T^ . The effects of the momentum spread of 
the distribution of the impurities on the collision rate is 
also examined by both analytical and numerical means. 

The paper is organized as follows. In Sec. ^ we in- 
troduce the physical system under study, while Sec. IIIII 
discusses the collision rate and reports analytical results 
in two limiting cases, the details being presented in an 
Appendix. In Sec. If VI we carry out Vlasov-Landau nu- 
merical simulations to examine the dynamics of the mix- 
ture. Finally, Sec. IVl offers some concluding remarks. 



II. THE PHYSICAL MODEL 

The system that we study is a spin-polarized Fermi gas 
in a trap containing bosonic impurities, which are free to 
move inside their own trap and interact with the fermions 
both through a mean-field potential and through colli- 
sions. The temperature of the system is finite and well 
below the Fermi temperature Tp , but well above the crit- 
ical temperature for Bose-Einstein condensation of the 
impurities. 

The fermionic component (j = F) and the bosonic 
impurities {j — B) in harmonic confining potentials 

V^xt (i") ^I'S described by the one-body distribution func- 
tions /^"'•'(r, p, t). These obey the Vlasov-Landau kinetic 
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equations (VLE), 

djU)+P_.yjU)^y^U^j).y fU)^C[f^F)j(B)^ (1) 
rrij 

where the Hartree-Fock effective potential is U^^^ (r, t) = 

V^xt (i") + gn^^'^ {r , t) with j denoting the species different 
from j. Here we have set g — 27r?i^ a/m^ with a being the 
s-wave scattering length of a fermion-boson pair and rrir 



where = f-^Hr,p,t) and /|^' = /(^^ (r, p„ i), d% is 
the element of solid angle for the outgoing relative mo- 
mentum P3 — P4, V — |v — V2I is the relative velocity of 
the incoming particles, and a = Ana'^ is the scattering 
cross-section. The collision satisfies conservation of mo- 
mentum (p + P2 = P3 + P4) and energy {e + e2 = £3+^4), 
with Ej =p]/2mj + U'^^l 

Hereafter we shall focus on a specific mixture of exper- 
imental relevance, namely, a '^''K-^^Rb gas with strongly 
attractive scattering length a = —410 Bohr radii. At 
this bare scattering length the mixture is not driven to 
collapse for the values of the atomic densities that we 
consider in this work 1^]. We shall focus on the case of 
isotropic traps with trap frequencies ojp 7^ ws, chosen as 
the geometric average of the experimental ones i.e. 
LUF ^27T X 134s~^ and lub = 27t x 91. 2s"^ 

III. THE COLLISION RATE 

The physical observable that identifies the dynamical 
regime of the gas is the quantum collision rate Tq. In 



Here and v.; are the total momentum and the rela- 
tive velocity of the incoming particles, and 6i {Of) is the 
direction of the relative velocity of the incoming (outgo- 
ing) particles in the scattering plane. At equilibrium the 
Fermi and Bose distributions in Eq. ||2J) are given by 

/(^■)={exp[/?(ep')-/i(^"))]+a-' (4) 

where ^ = 1 or — 1 respectively for fermions and bosons, 
/i^-'^ are the chemical potentials, and j3 = l/fc^T. The 



their reduced mass, and rS^^ (r, t) is the spatial density 
given by integration of f^^^ (r, p, t) over momentum de- 
grees of freedom. Since we deal with low concentrations 
of impurities we have neglected impurity-impurity inter- 
actions. In addition, collisions between spin-polarized 
fermions are negligible at low temperature and thus the 
collision integral C in Eq. Q involves only collisions be- 
tween fermions and impurities. This is given by 



(2) 

I 

the collisionless regime the atoms collide only rarely and 
we have Tq <^ cuj, while if F^ ^ ujj the gas is in a col- 
lisional regime which can be well described by hydrody- 
namic equations. In the general case a kinetic treatment 
is necessary such as the one used in this work. The colli- 
sion rate can then be evaluated either from a numerical 
simulation which actually counts the number of collisions 
at each time step, or by direct integration of the collision 
integral C over momenta. Numerical simulation runs for 
the problem at hand will be presented in Sec. IIVI 



In the second type of approach that we have mentioned 
above, the starting point is the usual assumption that be- 
fore a collision the distribution functions for the incoming 
particles and for the occupancy of the final states arc the 
equilibrium Fermi and Bose distributions. The local col- 
lision rate r|^°'^ is then given by 



(3) 

I 

total number of instantaneous collisions occurring in 
the system per unit time is then calculated as F, = 
/ F^°'^ d^r. 

One can further develop the evaluation of Eq. Q in 
two limiting cases (see Appendix). In the first case the 
temperature T of the fermions is very low - more pre- 
cisely, Tc T <$:Tf where Tp = {mFY'^fiujF /ke is 

1/3 

the Fermi temperature and — O.QANg hLOB/ks is the 
critical temperature for Bose-Einstein condensation, with 



I 

c = jdh2 dn^ V [(1 - /(-))(i + /f - /^"Vf ^(1 - fP)ii + f^], 



I 

rf = ^ I d^P. I dv.v! I dcose. I dcos^,/rVf ^(1 - fP)ii + ff') ■ 
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Np and Nb being the numbers of particles of the two 
species. In Eq. (jSJ we can then replace /'^^•' by the zero- 
temperature Fermi distribution and f'^-^^ by a classical 
distribution of the form 

fiB)(r,p) = ^l^ps{r)e-r''/Po (5) 

where PB(r) is the local density and pq the momentum 
spread of the impurities. 

The parameter po is crucial in determining the scat- 
tering rate in this situation. When the momentum 
spread of the impurities is low, i.e. for pq/pf <C 1 
where pp — \j2rap all collisions are forbidden 

by the Pauli principle since any fermionic final state 
which would be allowed by kinematics belongs to the 
Fermi sphere. Collisions involving empty fermionic fi- 
nal states become possible on increasing po and the col- 
lision rate starts increasing. At the same time the vol- 
ume occupied by the impurities in the trapped system 
also increases and reaches a maximum value when the 
two clouds have approximately the same size, i.e. for 
Po ~ "TfiB^BVFlijnpijjp). On an additional increase of 
Po the effective number of impurities that can interact 
with the fermions goes down and consequently the num- 
ber of collisions decreases. In this latter limit Pq/pp ^ 1, 
the leading term of the local collision rate is calculated 
in the Appendix and leads to the result 

T.c.^o'^^^NbNp. (6) 
7^ Po 

A more direct estimate in the same limit is obtained by 
a classical evaluation of the collision rate through the 
relation F^ ^ aviNpNs with the impurity density Nb ~ 
NB{mBi^B /po)^ / {4:n /3) . The relative velocity Vi can be 
set to Po /tob and we get 

r,^3a^^^NBNp, (7) 
Po 

which is in good agreement with Eq. ©. Equations © 
and ((TJ show that the collision rate rises for heavier impu- 
rities, thus favoring a coUisionless-to-coUisional dynami- 
cal transition. One may consider using in this context 
not only atomic impurities but also e.g. heavier stable 
molecules. 

The second limiting case in which progress can be made 
by analytical means refers to a gas at high temperature, 
it is worth noticing that in this case there would be no 
difference in considering as impurities either bosons or 
fermions. The behavior of the mixture approaches that 
of a classical two-component fluid, where the number Td 
of collisions per unit time is 

^3/2 

Tei = ^^NBNp(3 (8) 

(see Appendix). In Eq. (jSJ A:^ = niBi^B ^f^f/ (''^b^b^ 
mpujp) is the reduced oscillator strength. 



Finally, in close-to-equilibrium situations the evolution 
of the fermions can be obtained from a first-order approx- 
imation on the collision integral as 

C = -7(/(^)-/r) (9) 

f p\ 

where /q is the equilibrium distribution. The quan- 
tity 7 determines the damping rate of oscillatory mo- 
tions of the fermions and at low coUisionality an order- 
of-magnitude estimate of the damping can be obtained 
by evaluating one of the terms entering the full collision 
integral (see e.g. 0). This yields a relationship between 
collision rate and damping rate as 

j^r,/Np. (10) 

The accuracy of this relation will be explored numerically 
in the next Section. 



IV. DYNAMICS OF THE FERMIONS 

The analysis of the dynamics of the mixture requires 
the solution of the VLE in Eq. QJ. The solution is ob- 
tained by using a fully three-dimensional concurrent code 
with numerical procedures which are essentially the same 
as those described in Refs. 0, for a two-component 
Fermi gas, except for the choice of the initial phase-space 
distributions and for the bosonic enhancement factors en- 
tering the collision integral. The details have already 
been given elsewhere [isj . 

We shall analyze two different experiments. In the 
first, the mixture is prepared in an equilibrium state 
with distribution functions as given in Eq. In the 
second experiment the fermionic and bosonic clouds are 
prepared independently at two different temperatures 
(T and Tb) and are then superposed in space. In this 
case the fermions obey a phase-space distribution given 
by Eq. with a vanishing mean-field term, while 

the impurities are distributed according to Eq. (O with 
Po = (2mBrB)i/2. 

A. Dynamics near to equilibrium 

Having prepared the mixture in an equilibrium state 
inside isotropic traps, the system is set off by giving a 
slight displacement to the center of the fermion trap. 
We count the collisions made per unit time and per 
fermion, and report its average during roughly 8 oscil- 
lations. Therefore, at variance from what was discussed 
in Section lllll the collision rate calculated in this way 
includes both the effects of the instantaneous collisions 
and of the dynamics. 

In Fig. ^we report the scaled collision rate Tq/(3 and 
(in the inset) the collision rate F^ as functions of Nb /Np 
at Np = 10^, for three values of the temperature below 
the Fermi temperature. We are dealing in these cases 
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FIG. 1: Scaled collision rate Tq/PNp per fermion (in units 
of huip) in a ^"K-^^^Rb mixture in isotropic confinement with 
Nf = 10^, as a function of the impurity concentration Nb/Nf 
at three values of the temperature T/Tf (symbols). The 
dashed line shows the classical limit of Eq. The inset 

shows the collision rate Vq/NF per fermion, in units of ujf- 
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FIG. 2: Damping rate 7 (in units of u>f) in a «K-»^Rb mix- 
ture in isotropic confinement with Nf ~ W, as a function 
of the coUisionality Fq/NFi^F per fermion (cf. Fig. and 
at three values of T (symbols). The dashed line shows the 
first-order approximation given in Eq. IIUH . 



with almost classical impurities (T/Tc ~ 1.2 — 16) giv- 
ing a very small contribution to the mean-field potential. 
The coUisionality of the gas obviously decreases with Nb 
and goes linearly to zero for vanishing Nb ■ The effect of 
cooling the gas is instead more subtle. In spite of Pauli 
blocking the collision rate increases on cooling (see inset), 
as a result of the increase of the density of impurities in 
the overlap region of the two clouds. This effect has the 
same origin as the factor /3 entering the classical collision 
rate in Eq. JSJl and indeed the scaled quantity Tq/P de- 
creases on cooling, as is shown in the main body of Fig. 
n This also shows that the classical formula in Eq. ^ 
gives a good account of our results at T — 0.6 Tp (see 
also Ref. p3|)- Similar trends are obtained in simulation 
runs at A^F = lO'^. 

We have also analyzed the damping rate 7 of the oscil- 
lations of the fermion cloud under the influence of the im- 
purities and its relationship with the collision rate. With 
the present system parameters the gas is close to the col- 
lisionless regime, as is demonstrated by displaying in Fig. 
|21the relationship between 7 and Tq for isotropic confine- 
ment. The damping rate was obtained from the simula- 
tion data by fitting the center-of-mass coordinate of the 
fermions to the expression zp{t) cx cos{ujp t) exp{—^ t). 
It is seen from Fig. |21 that the first-order estimate given 
in Eq. (|10|l for a fiuid at low coUisionality is not far from 
our numerical results. For what concerns the dependence 
of 7 on temperature, we see from Fig.|2that it is closely 
related to that of since in this type of plot the data 
for different temperatures are essentially superposed. 

B. Dynamics far from equilibrium 

We focus on an experiment in which, instead of varying 
the number of impurities or the shape of the trap, we can 
control the distribution of the impurities. In practice this 



could be realized by preparing spatially separated clouds 
and then rapidly transferring them into the same spatial 
region. This situation is closer to the system discussed 
m Ref. 0, where strongly out-of-equilibrium collisions 
were assumed. Being an out-of-equilibrium situation we 
choose to characterize the dynamics by displaying the 
time-averaged collision rate of the system as function of 
Po- 

The collision rate per fermion as a function of the 
momentum spread of the impurities is shown in Fig. |31 
Here the time-averaged results from simulation runs on 
a Fermi gas at T = 0.3 Tp are compared with those from 
the numerical integration of Eq. ^ for a fully degenerate 
Fermi gas (T = 0) at equilibrium and from its analytical 
expansion at large po as given in Eq. I^. Even though 
there are quantitative differences between the results at 
T = and T ~ 0.3 Tp, the qualitative behavior is the 
same at both temperatures and the discrepancies could 
be attributed to the enhancement of the coUisionality at 
higher T and to out-of-equilibrium effects missing in Eq. 
P|l. The peak near pp corresponds to the situation 
where the two clouds essentially occupy the same spa- 
tial region. Correspondingly, in view of Eq. H10|l also the 
damping rate will have a peak at the same point. On 
further increasing the momentum spread po the effec- 
tive number of impurities that interact with the fermions 
diminishes and consequently the collision and damping 
rates decrease. 

Figure 0] shows the time-averaged ratio (Tq/Tci) be- 
tween the quantum and the classical collision rates as a 
function of Tb/Tp, as estimated during the simulation 
in correspondence to the data points in Fig. |21 Almost 
all collisions allowed by kinematics for po > pp involve 
empty fermionic final states so that Pauli blocking is in- 
effective. For pq < Pp instead, only impurities in the tail 
of the momentum distribution can contribute to the col- 
lision rate and the Pauli suppression increases dramat- 
ically. The quantum suppression is still lower than in 
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FIG. 3: Collision rate per fermion (in units of luf) as a func- 
tion of the momentum spread po (in units of the Fermi mo- 
mentum pf) for Nf = lO"* ""K atoms and Nb = 100 *^Rb 
impurities inside spherical traps. 
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FIG. 4: Time-averaged ratio (Tq/rd) between the quantum 
and the classical collision rate at T = 0.3 Tf, obtained in the 
numerical simulation as a function oi Tb/Tf (in log scale) for 
the gas shown in Fig. |3| 



vents these limitations and may offer a feasible method 
to increase the coUisionality of a spin-polarized fermion 
gas. 

For near-to-equilibrium dynamics in ^'^K-^^'Rb mix- 
tures similar to those used in actual experiments at LENS 
[T^ we have found that collisions are rare and that the 
gas is close to the collisionless regime. We have suggested 
that some parameters characterizing the system could be 
tuned to increase its coUisionality and thus drive the sys- 
tem towards the collisional regime. We have also shown 
how in a far-from-equilibrium experiment on a fermion 
gas at T ~ the addition of bosonic impurities with mo- 
mentum spread around the Fermi momentum pp could 
induce an enhancement of coUisionality. Due to the de- 
pendence of the collision rate on the mass of the impuri- 
ties, a further increase of coUisionality could be achieved 
by choosing as imp urities heavy particles such as ^^'^Cs 
|2Cl| and ^^^Yb [2l| atoms, or strongly bound molecules. 

Finally, the trap anisotropy could be raised so that 
for cigar-shaped traps the axial oscillations would be 
damped more rapidly than the radial ones on their own 
time scale. The system could thus be driven into the in- 
termediate collision regime on increasing the anisotropy 
of the trap. The transition from the collisionless to the 
collisional regime as driven by anisotropy will be investi- 
gated in detail in future work [2^ . 
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the case of a two-component mixture of fermions, where 
both components are subject to Pauli blocking. In this 
case (Tn/Tri) has been estimated to be roughly 0.65 at 
T = 0.3 Tp [13, while in the present case of bosonic im- 
purities we find (L,/Fc/) ~ 0.9 for Tb = 0.3 Tf (see Fig. 

EI) 

V. SUMMARY AND CONCLUDING REMARKS 

We have studied the collisionless properties of a spin- 
polarized Fermi gas interacting with a small number of 
thermal bosonic impurities and have focused on two re- 
lated aspects of its coUisionality, i.e. the collision rate 
and the damping rate of oscillations. While the superflu- 
idity of Bose-condensed atoms lowers the coUisionality of 
Bose-Fermi mixtures and the Pauli principle operating in 
a second fermionic component blocks collisional events in 
Fermi- Fermi mixtures, the use of thermal bosons circum- 



APPENDIX A: SOME ANALYTICAL RESULTS 

1. High temperature 

At high temperature the local collision rate in a two- 
component mixture with only inter-species interactions 
is given by the classical expression 

r':r{r)^a J dVd'p2^;/(^Hr,pi)/(2)(r,p2) (Al) 

where v is the relative velocity of particles 1 and 2 and 
the classical distribution functions in an anisotropic trap 
are 

/^^■Hr,p).^/?3^|e-K^--'^'«). (A2) 
Integration of Eq. ljAl|l is straightforward and gives 
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^ cl 



(r) 



2V2 7r7/2 



1/2 

rrir 



(A3) 



if mean-field potentials can be neglected. The total col- 
lision rate can then be calculated by evaluating the inte- 
gral of Eq. HA3|) over space. This yields Eq. © for the 
classical collision rate. 



2. Low temperature and large momentum spread 
of bosonic impurities 

The collision rate can also be analytically estimated 
when the fermions are at T ~ and for large values 



of the momentum spread po- In this case the fermionic 
distribution is the Fermi-Dirac step function, while the 
impurities are taken thermally distributed according to 
Eq. (0). Neglecting the Bose enhancement factor and 
mean-field effects, rj^°'^ can be written as 



T^IOC 
? 



(r) 



^ mim2 PB{r) 



mi 



m2Vif 



[12 



"0 e 



(A4) 



where the integration bounds in Vi, i.e. v- = \Pi/m2 — 
pp/mr\ and w+ — {Pi/m2+PF/'mr)^ come from requiring 
that the final states in a collision are unoccupied. 

To obtain the number of collisions per unit time we 
integrate the above expression inside the Thomas-Fermi 
radius i?TF of the fermion cloud. This is defined as the 
classical turning point of a fermion with energy /.t*^^^ in 
the external trapping potential, i.e. 

Rtf = \ 

y rapojF 

where we set /i*-^-* at the non-interacting value /i*-^-* = 
p\/2mp = {<o NpY^'^ hup. On increasing po the inte- 



grand in Eq. (jA4(l is well approximated by an inverted 
parabola with zeroes at u_ and , and hence the integral 
over Vi is trivial. We obtain 



^ 3 m2 TT^/^h 



Expression ^ in the main text follows from taking 
a local-density approximation for pi? by setting pp — 

\J 2mp{^^^'> — V^^^{r)) and performing the space inte- 
gration of the local collision rate in Eq. ljA6|l . 
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